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Abstract  - Tle  thermal  self -focusing  of  a Gaussian  electromagnetic  beam 
incident  upon  a dense,  low  temperature  plasma  is  discussed.  The  behavior 
of  the  beam  as  it  propagates  through  the  plasma  is  analyzed.  A measure 
of  the  penetration  depth  of  the  beam  in  an  inhomogeneous  plasma  as  a func- 
tion of  its  power  is  given. 


INTRODUCTION 


The  nonuniform  intensity  of  the  electric  field  in  a Gaussian  electro- 


magnetic beam  incident  upon  a plasma  can  cause  a redistribution  of  the 
electrons  in  the  plasma.  This  mechanism  leads  to  the  phenomenon  of 
self -focusing  and  enables  the  beam  to  penetrate  into  an  otherwise  over- 
dense  plasma  (Akhmanov,  Sukhorukov  and  Khoklov,  1968).  It  is  well- 
known  from  geometric  optics  that  a ray  incident  upon  a region  of  increasing 
index  of  refraction  will  refract  into  the  higher  region.  In  a plasma  the 
index  of  refraction  is  approximately  given  by 


P ‘ (wpe/w2)  * i(vei“pe/w3)ll/t"  = [(1  ' (ne/nec)  " ^eiV^ec^  2 
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where  n is  the  critical  electron  density,  i.  e.  the  plasma  density  where 
ec 

the  plasma  frequency  equals  the  incident  laser  frequency.  Thus,  if  a ray 
is  incident  on  a plasma  with  a density  gradient,  the  ray  will  refract  into  a 
less  dense  region.  Hence,  if  somehow  the  incident  electromagnetic  radia- 
tion can  induce  a change  in  the  initially  uniform  plasma  density  distribution 
in  the  transverse  direction,  then  self -focusing  may  take  place. 

In  the  case  of  a fast  (laser  pulse  time  shorter  than  the  energy  relaxa- 
tion time  of  the  electrons)  pulsed  laser  interacting  with  a strongly  ionized 
collisionless  plasma,  the  ponderomotive  force  is  able  to  redistribute  the 
electrons  in  such  a way  that  there  exists  a positive  density  gradient  in  the 
transverse  direction  away  from  the  beam  axis.  Thus  a transverse  gradient 
of  the  index  of  refraction  is  established,  and  this  leads  to  self -focusing  of 
the  incident  laser  beam. 

However,  in  the  case  of  a low  temperature,  high  density  plasma  such 

as  that  employed  in  our  laboratory  Z-pinch  experiment  (Rockett,  et  al. , 

12 

1978),  the  electron-ion  collision  frequency  is  very  high,  v . = 3 x 10  sec. 

ei 
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Moreover,  the  pulse  time  of  our  COg  laser  is  of  the  order  of  3.5  x 10  sec. 
i.  e.,  much  longer  than  the  energy  relaxation  time  of  the  electrons.  There- 
fore, in  this  case,  the  nonlinearity  appearing  through  the  nonuniform  heating 
of  the  electrons  by  the  incident  laser  beam  is  more  important  than  that  due 
to  the  ponderomotive  force  effect,  and  we  have  a thermal  self -focusing 
phenomenon.  This  is  characterized  by  a redistribution  of  the  plasma  elec- 
trons by  the  nonuniform  heating  in  a direction  normal  to  the  incident  beam 
axis.  The  transverse  nonuniform  distribution  of  the  electric  field  intensity 


of  the  incident  laser  beam  creates  a nonhomogeneous  electron  temperature 
in  the  radial  direction  away  from  the  beam  axis.  For  a Gaussian  electric 
field  in  the  laser  the  electrons  on  the  beam  axis  become  hotter  than  those 
away  from  it.  This  also  can  cause  self -focusing  of  the  incident  beam  as 
will  be  discussed  below.  (Sodha,  Ghatak  and  Tripathi,  1976). 


THERMAL  SELF-FOCUSING 

The  energy  of  the  electrons  acquired  from  the  incident  laser  beam  is 
dissipated  through  collisions  with  the  ions  and  through  thermal  conduction. 
The  rate  of  energy  lost  due  to  electron -ion  collisions  may  be  described  by 


3T<fkBVe-i  = !V2m'eq(VTi>/M> 
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where  v is  the  equilibration  time  of  the  electron  and  ion  temperatures. 

On  the  other  hand,  the  rate  of  energy  loss  per  electron  due  to  thermal  con- 
duction is 


dt  (2  kB  Te*tc 


1 9,  8Te, 
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where  k is  the  coefficient  of  electron  thermal  conduction. 

For  our  laboratory  Zrpinch  plasma  experiment  the  ambient  conduc- 
tions are  described  by 
T 


eo 


= 2.  32  x 105  °K. 


T.  < T 
10  ~ eo 

g 

Ve  = electron  thermal  velocity  <=»  4.  2 x 10  cm/  sec. 

in10  -1 
v =*  10  sec 

eq 
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v = electron-ion  collision  frequency  =>  3 x 10  sec 
ei 


r = radius  of  focal  spot  of  the  incident  laser  beam  «=  6 x 10  cm. 
o 

For  these  experimental  values  it  may  be  shown  that  the  energy  loss  of  the 
electrons  due  to  thermal  conduction  is  much  larger  than  that  due  to  elec- 
tron-ion collisions.  An  estimate  of  the  relative  importance  of  the  two  dis* 
sipation  mechanisms  yields 

(3/2)k  2mv  (T  -T.)/M  i„  v v . , 

1__B «■)'  « jL  . 3 22.  ° *2  Si  a o(10 -2) 

B , lie,,  M v2  ' 

Tt  (r  « ~)/ner  e 

for  our  laboratory  experiment.  This  suggests  that  we  have  a thermal  fo- 
cusing process  in  which  the  main  energy  dissipation  mechanism  is  elec- 
tron heat  conduction  in  the  direction  transverse  to  the  axis  of  the  incident 
beam. 

Next  we  calculate  the  electron  density  distribution  in  the  transverse 
direction.  The  steady  state  energy  balance  equation  may  be  written  as 

Id  d T 

e E • V = (nfi  rf  ^(r^)  (1) 

Substituting  the  velocity  from  the  equation  of  motion  for  the  electrons  un- 
der the  influence  of  the  oscillating  electric  field  of  the  CO^  laser  beam, 
we  obtain 


n v . e EE*  . , d T 

e ei  1 d e . 

2 “ " r dr  <rK  'dr- * (2) 

m w 

where  w is  the  frequency  of  the  incident  CC^  laser  beam,  and  we  have  used 
2 2 

the  fact  that  v . <<  u>  in  our  experiment.  We  now  assume  a Gaussian 
ei 

distribution  for  the  incident  laser  beam: 


EE*  = KZq  exp[-  r2/r¥(z)]/f2(z) 
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where  Eq  is  the  amplitude  of  the  electric  field  of  the  incident  beam,  and  f(z)  is 

a variable  which  parametrizes  the  beam  width. 

In  the  linear  approximation  we  carry  out  a Taylor  series  expansion  of  the 

right  hand  side  of  eq.  (3)  with  respect  to  r = 0,  i.  e.,  the  beam  axis.  Also  we 

use  ambient  values  for  the  coefficients  in  the  energy  balance  equation  (2).  We 

then  solve  for  T (r),  and  obtain 
e 


Te,r|  = Teo+  <Mo/4*of2|(rf[1-<I/«2U  - r2(l  - (r24r¥)]} 


, „ _ 2,2 

where  p = n e v . /mu 
o eo  eio' 

and  the  subscript  zero  denotes  ambient  values. 

The  sound  speed  in  the  z-pinch  plasma  in  our  experiment  is  about  2.  2 x 
6 . 

10  cm/sec.  As  the  characteristic  length  of  the  laser  beam  is  of  the  order  of 


(4) 

(5) 
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6 x 10  cm  (i.e.,  the  focal  spot  radius),  the  characteristic  time  scale  t 


H 


for  hydrodynamic  motions  transverse  to  the  laser  beam  axis  is  of  the  order  of 


, -9 

3 x 10  sec.  As  indicated  previously  the  laser  pulse  time  is  approximately 

-8 

3.5  x 10  sec.  Thus,  in  the  time  scale  of  the  laser  pulse,  i.e.,  for  t > t , 

H 

the  plasma  pressure  pg  + p.  is  approximately  constant  across  the  laser  beam. 

Furthermore,  ng  a n.,  and  a T.  as  the  time  scale  for  the  electron  and  ion 

temperatures  to  equilibrate  is  much  less  than  t . Consequently,  p = n k T 

H e e B e 

=s  constant  across  the  beam.  This  implies  that 

n (r)  sr  n T /T  (r)  /6) 

e eo  eo2  e 

where  ng  = n^  and  Tg  = T^,  the  ambient  values,  when  r = r^  From  T (r) 
as  expressed  in  eq.  (4)  we  then  obtain  the  electron  density  distribution: 


ne(r)  = "eo  - neo[{(ro  ’ ^ ^ } /{■ (4k  0Teof^  0Eo] 1 

+ [(rf  - r2)  - (l/4r2f2)(r*  - r4)]} ] ; 0 < r < r 


(7) 


This  expression  for  the  density  distribution  is  now  substituted  into  the  formula 
for  the  dielectric  constant  (neglect  the  imaginary  part  which  accounts  for  the 
absorption) 

e(r)  = 1 - [u>p(r)/u>2]  (8) 

The  result  is 


e(r)  = 6o  + (W2o/u,2)({(r2  -r2)  - (l/4r2f2)(r4  - r4)}/{(4*  ^f2/^2) 
+ [(r2  - r2)  - (l/4r2f2)(r^  - r4)]  }) 


(9) 


where  e = 1 - (w  /u  ),  i.  e. , the  dielectic  constant  of  the  ambient  plasma, 
o po  " 

Equation  (9)  gives  the  radial  variation  of  the  dielectric  constant  necessary 
for  the  thermal  focusing  process.  It  may  also  be  written  in  the  form 


e(r)  = eo+4>(Eo) 


(10) 


where  4>(Eq)  is  given  by  the  second  term  of  the  right  hand  side  of  eq.  (9).  It 
is  the  modification  of  the  plasma  dielectric  constant  by  the  electric  field  of  the 
incident  beam  through  the  nonuniform  heating  of  the  electrons. 

In  the  last  section  we  solved  the  plasma  heat  balance  equation  by  ex- 
panding the  Gaussian  beam  intensity  profile  about  its  centerline.  This  approx - 

2,22 

imation  is  strictly  not  valid  when  r r f is  of  order  unity.  Therefore  we 

o 

solve  the  heat  balance  equation  numerically  by  finite  difference  method  as 


follows. 


7 


First  we  note  that  the  collision  frequency  and  electron  thermal  conduc- 
tivity vary  with  the  temperature  as  (Spitzer,  1956) 


v . = v . [T  /T  (r)]' 
ei  eioL  eo  e 


k = /c  [T  (r)/T  ]" 

e eoL  e ' eoJ 


where  we  have  assumed  that  the  electron  pressure  is  nearly  constant  across 


the  laser  beam  (see  the  previous  discussion).  The  subscript  zero  indicates 
ambient  conditions.  The  heat  balance  equation  may  now  be  written  as 


(P„TLA_TI/2)EE*(r)  = -“(r  Tf/2aTyar) 


o eo  eo  e 


2 2 

where  = nQ  vei0e  /mu  * Using  the  Gaussian  beam  intensity  profile  and 

7/2 

defining  a new  variable,  F = T , we  can  rewrite  this  as 

(p  T6  EZ/kl  f2)F(r)exp(-  r2/r2f2)  = -|---^-(r^) 

' ~ o eo  o'  eo  r 'o'  7 r dr  dr  ' 

By  letting  R = r/r^f,  this  becomes 

2 1 dF  d2F 

- X F(R)exp(-  R ) = | (I! 

dr 

7 6 2 2, 

where  X = — P QT  E^r^/jc  eo’  Equation  (7)  is  subject  to  the  boundary  condition 

F(r  ) = T7^2  (12 

o eo 

We  solve  eq.  (11)  numerically  using  a finite  difference  scheme.  The  R varia- 
tion is  divided  into  N intervals  of  length  aR;  this  results  in  the  following 
finite  difference  equation  for  F: 


Fn+1  = {l/2n+l)}{4nFn  - n(AR)2XF^exp(-n2AR2)  - (2n-l)Fn_1}  (13) 


First  we  estimate  Fq*  Then  Fj  is  obtained  using  the  heat  balance  equation 

above  (valid  for  iR  <<  1).  With  these  two  starting  values  for  F and  F , 

o 1 

the  succeeding  F^  (n  > 1)  is  generated  using  eq.  (13)  above.  The  value  of 
F^  so  obtained  is  then  compared  with  the  boundary  condition  (12),  and  the 
next  estimate  for  Fq  is  adjusted  accordingly.  An  iteration  process  is  now 
employed  until  F^  is  as  close  to  T^  as  desired.  The  electron  tempera- 
ture distribution  in  the  radial  direction  obtained  by  this  numerical  scheme 
is  given  in  figures  1 and  2.  As  a comparison,  the  corresponding  electron 
temperature  distribution  calculated  by  expanding  the  Gaussian  incident  beam 
profile  is  also  shown.  Note  that  the  linearized  method  yields  unphysical  re- 
sults for  the  electron  temperature  distribution  as  shown  in  figure  2.  However, 
for  values  of  the  incident  beam  electric  field  intensity  not  far  above  threshold, 
the  linearized  method  yields  reasonable  results. 

Another  way  of  solving  the  energy  balance  equation  (2)  is  given  in  the 
Appendix.  This  procedure  is  an  improvement  upon  the  linearization  method. 

In  fact,  it  is  an  e".-»ct  method  of  solution  except  for  the  initial  integration  of 

the  energy  balance  equation,  where  the  coefficient  8 = ne^v  ./ mw  2 is  assumed 

ei 

to  be  independent  of  the  electron  temperature.  As  seen  in  the  Appendix  this 

method  also  leads  to  an  expression  similar  to  eq.  (9)  for  <t>(E‘’). 

o 

The  fundamental  equation  governing  the  propagation  of  an  electromagnetic 
wave  is 

- V[V-  E)  = w2  e E/  C2  (13) 

The  second  term  on  the  left  hand  side  is  negligible  in  comparison  to  the  first 

term  if  k ■v 


(In  e)  <<  1,  where  k is  the  wave  vector.  The  above  inequality 


d J I 


innppmpmn 


TF"”?1 
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is  satisfied  in  almost  all  cases  of  practical  interest  (Sodha,  Ghatak  and 
Tripathi,  1976).  For  a cylindrically  symmetric  electromagnetic  beam  we 
then  have 


2 

8 "E 


8 E 


8 z 


1 8 v ~ 2 2 

2+78l(r87)  = w e £/C 


(14) 


To  solve  this  equation  we  follow  the  method  employed  by  Akhmanov,  Sukhoru- 
kov,  and  Khokhlov  (1968)  using  the  paraxial  ray  approximation  in  geometric 
optics  (see  also  Sodha,  Ghatak  and  Tripathi,  1976).  We  first  write 


€ (r,  z)  = e (r=0,  z)  + 4>  (<  E > ) 


(15) 


where  ®(<E  > ) is  the  dielectric  response  to  the  incident  electromagnetic  beam 

2 

given  by  eq.  (12).  We  assume  that  <M<E  > ) <<  e(r  = 0).  Using  a WKB  approx- 
imation we  may  express  gin  the  form: 


E = A(r,  z)  [ e:1'  4(0)/c1/4(z)]  exp{i(wt  - f dz^  e l/  " c)} 


(16) 


with 


A(r,z)  = AQ(r,z)exp[  • i S(r,  z)]  (17) 

and  S is  the  eikonal  of  the  wave. 

Substituting  eqs.  (15)  and  (16)  into  eq.  (13),  making  use  of  eq.  (14),  and 
separating  the  real  and  imaginary  parts,  the  following  solution  for  a Gaussian 
incident  beam  is  obtained  (Sodha  and  Tripathi,  1977): 


S = j Y(z)r2  + ^(z) 

<5  j ( z)  = an  arbitrary  function  of  z to  be  determined  by  the  boundary  conditions 

2 2 2 2 2 2 

Ao  = <Vf  )exp(-  r >V  » (IS) 

1 / 2 , 1 

y = (u>  c /c)  — df/dz 


where  the  beam  width  parameter  f is  governed  by 


2 

d“f 

-> 

dr 


1 K df 


,3  £ 


(19) 


£ = z ''k  r is  a dimensionaless  "distance" 

o 

2 1/2 

a = (rQ  e f /2w)de/dz 

_ 2 2 4 2 

K - uj  E Sr  /4c  * T is  the  "self -focusing  constant." 

po  o o eo 

The  initial  conditions  on  f at  £ =0  are  f=l  and  df  ,'d|  =0.  This 
corresponds  to  an  incident  plane  wave  front. 

(a)  Homogeneous  plasma 

Let  us  first  consider  propagation  in  a homogeneous  plasma;  in  this 

case  the  last  term  on  the  right  hand  side  of  eq.  (19)  vanishes.  For  a 

2 2 

plane  wave  incident  at  £ =0  we  have  d f/d£  =■  0 everywhere,  if  the  first 

and  second  terms  on  the  right  hand  side  of  eq.  (19)  cancel  each  other, 
df 

Thus,  — will  remain  at  zero,  and  hence  f=l  for  all  values  of  £ . In  this 

d§ 

case  the  beam  propagates  without  modification,  neither  diffracting  nor 
focusing.  This  condition  thus  yields  the  threshold  intensity  for  thermal 
self -focusing: 


E = 4c^k  T /w  l3  r^ 

o crit  B eo  po H o 


(20) 


The  behavior  of  the  beam  as  it  propagates  through  the  plasma  (2)  may  be 
obtained  by  solving  eq.  (19).  A numerical  analysis  of  eq.  (19)  yields  the 
variation  of  the  beam  width  parameter  f with  the  dimensionless  "distance' 


£ . The  result  is  shown  in  fig.  3. 


r 


•i 


* 
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The  self -focusing  constant  K can  be  related  to  the  incident  beam  power 

by  integrating  the  expression  for  the  Poynting  vector 

, oo 

/01  i „ .2,2. 


P = 


= (eo  c/8t r)J  dr  2Trr(EVr)exp(.  ^/rV) 
o ° 

}/2 

o 

Hence,  we  may  write 


ej  c EZrZ/8 


o o' 


K = 2u5po-.3ro  P/°3  *1qZ  K T 


eo 


(21) 


Thus  K is  directly  proportional  to  the  incident  laser  power  P. 

The  critical  power  for  thermal  self -focusing  i6  obtained  from  eq.  (20) 


as 


p rit  = c3e1/Z  kT  /2co2  0r2 


eo7  pcA  o 

Using  Spitzer*  s values  for  vg.  and  K this  can  be  approximately  expressed 


as 


39  _5  ,3 


Pcrit  - 4*  3 x 10  T*/n„  ergs/sec. 


e'  e 


With  the  laser  power  of  3.  7 x 1014  ergs/sec  used  in  our  laboratory  experi- 


ment, the  critical  density  is  calculated  to  be  n = 2 x lO1^  cm-2.  Thic 


ec 


is  well  satisfied  by  our  Z-pinch  plasma. 

A parameter  of  interest  is  f^,  the  minimum  value  attained  by  f as 
the  beam  propagates  through  the  homogeneous  plasma.  This  parameter 
signifies  the  extent  of  focusing  associated  with  a particular  value  of  K. 

It  is  found  by  integrating  eq.  (19)  with  « = 0 once  to  obtain 


(5|)2  = 1 _ (1/^)  - 2K  In  f 


(22) 


P 
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I 
1 

The  minimum  value  of  occurs  when  df/d£  = 0.  Figure  4 shows  the 
behavior  of  f^  as  a function  of  K. 

(b)  Inhomogeneous  plasma. 

For  the  case  when  the  plasma  has  a positive  density  gradient  in  the 
direction  of  propagation  of  the  incident  laser,  we  can  calculate  the  penetra- 
tion distance  of  beam.  For  simplicity  we  consider  the  case  of  a linearly 
increasing  density  profile  of  the  form 

%„<*>  1 %[l  + (s/L>J  (23) 

where  (i  ^ is  the  electron  density  at  the  position  where  the  incident  laser 
impinges  upon  the  plasma,  i.  e.,  the  edge  of  the  cylindrical  plasma  column 
in  our  laboratory  experiment.  Lisa  density  gradient  scale  length. 

As  the  beam  penetrates  into  the  plasma  with  increasing  density,  the 
axial  dielectric  constant  decreases  and  a turning  point  is  expected  where 
e(r=0,  z)  vanishes.  However,  our  method  of  analysis  is  not  applicable  at 
this  point.  Therefore,  in  order  to  have  an  idea  of  the  penetration  distance 
of  an  incident  laser  beam  as  a function  of  the  power  of  the  beam,  we  intro- 
duce, albeit  arbitrarily,  a typical  penetration  distance,  z , where 

P 


e(r-0,  z^)  = 0.  69.  At  this  position,  the  ambient  plasma  density  at  r=r  is, 
of  course,  larger  than  the  density  corresponding  to  e = 0.69  at  the  beam 
axis.  The  variation  of  z^  with  the  incident  beam  power  is  shown  in  fig.  5. 


trie  field  strength  of  the  incident  laser  beam.  Thus  a beam  of  higher  power 

will  penetrate  deeper,  and  in  this  way  a sufficiently  intense  beam  can  penetrate 
into  an  overdense  plasma. 

Note  - Recently  Sharma  and  Tripathi  (1978)  made  an  attempt  to  interpret  the 
experimental  results  obtained  by  Rockett,  et  al.  in  our  laboratory.  However, 
they  assumed  an  incorrect  laser  beam  profile,  namely,  one  which  peaked  at 
r = rQ  > 0.  Our  observations  clearly  indicates  that  the  focussed  laser  beam 
has  a maximum  at  r = 0,  i.  e.  , the  axis  of  the  laser  beam.  Hence,  a Gaus- 


sian profile  is  the  more  correct  one  to  assume. 
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Appendix  - 


The  energy  balance  equation  (2)  may  be  written  as 


(7pE2/20f2)exp(-  r2/r¥)  = . 1 ^_{r  ^(T7/2)} 

° o r dr  dr  e 

where  the  coefficient  of  thermal  conduction  is  written  as  k = 6 T^/  2,  with 

0 a constant  (Spitzer,  1965).  Also  (3  3 nv^e  • We  now  assume  that 

P = P0*  where  the  subscript  zero  indicates  ambient  conditions.  Then  (Al) 

may  be  integrated  once  to  yield 

(cj/r)  + (pr2f2/2r)exp(-  r2/r2f2)  = ^-(T^/2) 
where  p = 7(3E2/20f2. 

As  dT  /dr  = 0 at  r = 0,  therefore  c = - pr2f2/2.  Hence 
c l o ' 

* (p  rff2/2r)  + ( P r 2f2/ 2r)exp(  - r2/rV)  = ^_(T7/2) 
u o o dr  e 

Before  we  integrate  (A2)  let  us  consider  the  temperature  dependence  of  the 
coefficient  (3.  As  p depends  on  the  temperature  through  the  density  and  the 
electron-ion  collision  frequency,  we  let 


= e0<VT.>7/2  * 

Equation  (A2)  may  then  be  written  as 


(Plrff2/r){exp(-  r2/r2f2)  - i}  = ~(T^) 
This  equation  can  be  integrated  formally.  The  result  is 


C2  +lPlrof2{ 


2 4 6 

r r r 7 

~TT  + “TT r5  + *“)  = TJr> 

rf  4r  f 18r  f e 

o O o 


.it 
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where  we  have  used  the  formula 


/(l/x)exp(ax)dx  = In  x + + §-~  + + . . . 

Applying  the  boundary  condition  T (r  = r ) = T we  obtain 

e o eo 

T«<r)  = Teo(,+<^„E>%T  r,((rf.r2,  -(l/4r¥),r4-r4  + 


+ (l/18rV)(/-r6)  - ...]  }1/7 


(A3) 


Proceeding  as  before  this  leads  to 


4>(Eo)  = -P/U  + (7P0Eo/4,coTeof2)l{ro"r2) 


- (l/4rV)(r4-r4)  + ...]  } 1/1 )) 


(A4) 


If  we  carry  out  a binomial  expansion  of  the  denominator,  we  obtain  the  pre- 

■> 

vious  expression  for  4>(E^)  as  shown  in  eq.  (9). 


^ .1 


Figure  Captions 


17  3 

Electron  temperature  versus  beam  radius,  n » 2 x 10  /cm  ; 

eo 

T ■ 20  eV;  E“  ■ 2.66  x 10^  ergs/cm^;  Incident  beam  is  at 
eo  o 

critical  power  for  thermal  self-focusing.  The  solid  curve 
corresponds  to  results  from  linear  theory,  and  the  dotted  curve 
corresponds  to  the  numerical  results. 

2 9 3 

Electron  temperature  versus  beam  radius.  Eq  - 4.26  x 10  ergs/cm  , 

i.e.,  beam  above  critical  power,  n and  T same  as  in  fie.  1. 

eo  eo  6 

Note  that  the  linear  analysis  yields  an  unphysical  temperature 

distribution  near  r/r  - 1. 

o 

Variation  of  the  beam  width  parameter  f with  distance  into  plasma 

for  different  values  of  the  self-focusing  strength 
4 2 2 2 

K ■ Br  w R /4C  <T  . K » 1.0  corresponds  to  the  critical  power, 
o p o eo  r 

The  minimum  beam  width  parameter  f versus  the  self-focusing 

m 

strength  K. 

Penetration  distance  versus  beam  intensity,  where  a ■ wp  6 r /4C  kT 

o o o eo 

17  3 2 

nor.  “ 2 x 10  /cm  ; T ■ 20  eV.  aE  is  the  self-focusing  strength  K 

cu  CO  O 


of  the  laser  beam. 


(r)-T.)/X 


Ka0.5 


